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Abstract 

Thermal fluctuation of the cantilever position sets a fundamental limit for the precision of any 
Scanning Force Microscope. In the present work we analyse how these fluctuations limit the deter- 
mination of the resonance frequency of the tip-sample system. The basic principles of frequency 
detection in Dynamic Scanning Force Microscopy are revised and the precise response of a typical 
frequency detection unit to thermal fluctuation of the cantilever is analysed in detail. A general 
relation for thermal frequency noise is found as a function of measurement bandwidth and can- 
tilever oscillation. For large oscillation amplitude and low bandwidth, this relation converges to 
the result known from the literature, while for low oscillation amplitude and large bandwidth we 
find that the thermal frequency noise is equal to the width of the resonance curve and therefore 
stays finite, contrary to what is predicted by the relation known so far. The results presented in 
this work fundamentally determine the ultimate limits of Dynamic Scanning Force Microscopy. 



I. INTRODUCTION 



Since its invention almost 20 years ago, Scanning Force Microscopy (SFM)[1J has become 
an extremely powerful tool for a huge variety of nanoscale investigations. With respect 
to resolution and sensitivity, Dynamic Scanning Force Microscopy (DSFM) [2] seems to be 
the most promising technique. Even though (true) atomic resolution was first achieved with 
contact mode in liquid environment j3], now atomic resolution studies are generally performed 
with DSFM working in UHV environment [IH6] . Recently even sub-atomic resolution has 
been reported using DSFM [7]. In spite of the impressing advances of SFM and DSFM 
we believe that the ultimate limit of these techniques is still an open issue. For most 
applications, the temperature at which quantum limits become relevant, Tq = huo/k ~ 
lfiK, is well below typical temperature ranges used in SFM. Correspondingly, either thermal 
vibration of the cantilever [BHTT] . or fundamental limits of the detection technique used[TTl- 
[13] - related to shot noise of the "detection particles" - determine the resolution in SFM 
and DSFM. In most practical applications, the fundamental limit of SFM and DSFM is set 
by thermal noise. Thermal noise in an SFM set up is the consequence of the equipartition 
theorem, which relates the mean energy of the cantilever with the thermal energy kT, 

\ c a th = \ kT (!) 
where aj h is the (mean square) displacement of the cantilever induced by thermal fluctuation, 
and c is its force constant. For high resolution distance measurements stiff cantilevers (HI US] 
should be used (5z = yJkT/c, with 5z fluctuation of tip-sample distance), while for high 
resolution force measurements soft cantilevers [TU] are needed (5F th = c ■ a t h = Vc kT). 
Although the mean displacement and the mean force fluctuation of the cantilever are im- 
portant quantities, in many applications they do not directly determine resolution, either 
because appropriate filtering significantly reduces the measured noise, or because the SFM 
technique used -as for example DSFM- is not directly limited by the displacement or the 
force. 

In typical DSFM applications the cantilever is excited at or near its natural frequency 
and the variation of its resonant properties -oscillation amplitude, resonance frequency or 
quality factor - are recorded. For imaging applications in liquids or air usually the oscillation 
amplitude is used as control parameter that defines constant tip-sample distance and the 
frequency shift is measured as secondary channel. In UHV applications DSFM operation is 

2 



the other way, the frequency shift is the control parameter for tip-sample distance and the 
oscillation amplitude -dissipation energy- is measured as secondary channel. At the moment 
DSFM is the most sensitive SFM technique to measure tip-sample interaction, which is 
detected as a shift of the resonance frequency of the system. Therefore, with regard to the 
ultimate limits of SFM a key issue -to be discussed in the present work- is to understand 
in detail the effect of thermal fluctuations on DSFM detection schemes and in particular on 
frequency detection. At present, the thermal noise density of a frequency measurement is 
assumed to be [10 HH I2D] 



with Q the quality factor, vq the resonance frequency, a the (root mean square!) oscillation 
amplitude and bw the bandwidth of the measurement. Note that this relation diverges 
for small oscillation amplitude and is proportional to the square root of the measurement 
bandwidth. 

In the present work we will revise in detail how thermal fluctuation limits the measure- 
ment of frequency shift. We find that neither for large bandwidth measurements nor for 
low oscillation amplitude relation [2] is correct. We present a general relation valid for all 
ranges of amplitude and bandwidth that agrees with relation [2] in the low bandwidth and 
large amplitude range. Finally, we present experimental data that unambiguously proves 
the validity of the general relation obtained in this work. 

II. FREQUENCY DETECTION IN DYNAMIC SCANNING FORCE MI- 
CROSCOPY 

A thorough discussion of frequency detection and DSFM is out of the scope of the present 
work (see the excellent original works [ 171 fT8] or recent reviews by Garcia and Perez[19j and 
Giessibl[20j), nevertheless we think it is important to revise some of its basic principles. 
The focus of this revision is rather on a comprehensive physical explanation of the technique 
than on a profound analysis of the underlying physiscs and statistical mechanics (see[18j) 
or the electronic details of its implementation(see [TSJ |2"T]). 

An externally driven SFM cantilever -usually driven by inertial forces- is a textbook 
example of a harmonic oscillator. In many aspects SFM and modern gravitational wave 




(2) 
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detectors are governed by the same basic principles (see, for example, [22]). As described in 
detail in appendix A, the dynamic properties of such a system are described by the quality 
factor Q, by a driving amplitude a$ and by its natural angular frequency u>o = 2ttu = 
y/c/m e ff, with vq natural frequency and m e ff effective mass of the cantilever |23j. In the 
steady state regime, the response of the system to a harmonic excitation a^UQ cos(27rz/£) 
can be described by a (complex) amplitude A{y), or by two components X{y) and Y(u), 
the in-phase and out of phase components of the oscillation, corresponding to the real 
and complex parts of the complex amplitude A{y). The time response of the system is 
then a(t) = X(y) cos(27rz/t) + Y (u) sin(27rz/i) = \A(u)\ cos (2irut + (p(u)) where <p(v) is the 
phase between the driving force and the response. At its natural frequency i/q the phase is 
— 7r/2, the (complex) oscillation amplitude is A(vq) = — iY(vq) = —ia^Q and the in-phase 
component X(vq) vanishes. 

Figure 1 shows a schematic set-up of the main components used in a DSFM detection 
electronics. The multiplication stages together with the filters essentially calculate the two 



components X{y) and Y(y) -relations 18 and 19 in appendix A- from the measured oscillation 
of the cantilever a(t). When enabled, the Pi-controller of the DSFM detection electronics 
adjusts the driving frequency vq of the excitation signal in order to have X(v) = 0. In this 
way, the system is locked to the natural frequency of the cantilever, tracks this frequency if it 
varies due to tip-sample interaction and generates an output proportional to the shift of the 
resonance frequency. In addition to the main DSFM components figure 1 shows the signals 
(in the frequency domain) along the different points in the DSFM detection path. A key 
component of any DSFM detection unit is the voltage (VCO) or numerically (NCO) con- 
trolled oscillator that generates the excitation and reference signals for the lock-in detection 
scheme. In most applications this oscillator drives the piezo that induces motion of the can- 
tilever. The corresponding deflection a(t) is analyzed by multiplication with two reference 
signals in quadrature. We will first assume the most general case where the frequency v of 
the cantilever motion and that of the reference oscillator are different. Then, multiplication 
of the deflection signal a(t) = ao cos(27rz/t) with the reference signal a re f(t) = cos(27ri/ re jt) 
results in two quadrature signals x q (t) = a(t) cos {2itv re ft) and y q (t) = a(t) sin {2nv re ft) with 
frequency components at v& = v — z/ re j, and = v + v re f- As shown in appendix A, x q (t) 
and y q (t) can be obtained from X(v) and Y(y) by multiplication with (Ma(£) + Ms(t)), 
where Ma(£) is a rotation matrix generating a clockwise rotation with frequency and 
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Ms(i) is a rotation matrix generating a counter-clockwise rotation with frequency z/£. In the 
frequency domain the signal a(t) is thus splitted and shifted to the sum and to the difference 
frequency (see figure 1). After the multiplication stages the two quadrature signals x q {t) 
and y q {t) are low-pass filtered over a timespan proportional to the time constant r of the 
filter. The precise time domain signals are determined in appendix A. Again, two matrices 



Mta(£) and Mrs(t) can be defined (relations 24 and 25 in appendix A), corresponding to 
a clockwise and counter-clockwise rotation with the frequencies and z/^. As compared to 
the first matrices, these "filtered" matrices have delay angles <^a ( r ) and <^s if) as well as 
multiplicative factors 1/(1 + An 2 r 2 v\) and 1/(1 + 47r 2 r 2 i/|). 

Usually, in DSFM the signal entering the DSFM detection unit is at the natural frequency 
uq, and the reference is at the same frequency v re f = vq. Then multiplication of the input 
signal with the reference signals will result in spectra around v = and v = 2v§. In this 



case Mta(^) = Id/2 , with Id the identity matrix, and we find, using relations 24 and 25 



. ... X(v) A(i/)(cos(2o;o^) — 2w t sin(2ci;oi)) F(^)(sin(2u;o£) + 2ojqt cos(2u;ot)) 
{x q {t)) T = — + 2(l + 4 W 2 r 2 ) + 2(l + 4 W 2 r 2 ) 

(3) 

Y(v) X(z/)(sin(2co> t) + 2u t cos(2u t)) Y(u)(— cos(2oj t) + 2u t cos(2u t)) 
{y q {t)h = — + 2(l + 4w 2 r 2 ) + 2(l + 4w 2 r 2 ) 

(4) 

These signals can be conveniently represented in the frequency domain (see fig. 1), where 
the first peak at v = is one-sided (since no negative frequency exists) and the second 
at v — 2z/ is two-sided. In typical DSFM applications the filter is usually set so that 
b>o » 1/r, then it blocks the 2i/o component and passes the signals from DC to Vq w 1/t. 
The value of r determines the speed and the "cleanness" of the signals. Large time constants 
(small bandwidth) results in clean but slow response, while, conversely, small-time constants 
will result in "unclean" signals - in particular with significant 2z/ component - but with 
fast response. In our system we have found that time constants of 3/uq to 10/b> give the 
optimum compromise between speed and "cleanness". 

In the frequency domain, the Fourier transforms of the signals x q (t) and y q {t) are simply 
multiplied by the frequency-dependent gain of the filter. Depending on the time constant of 
the filter, the total amount of signal may be decreased. The calculation of the total signal 
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Aubw (v c ) measured around a frequency u c within a certain bandwidth bw is performed most 
conveniently in frequency space: 



where v(y) is the (spectral) signal density (unit: V/ V Hz) and bw = 1/r the effective 

bandwidth of the filter. The signal densities corresponding to (x(t)) T and (y(t)) T in the 
frequency domain are X q {v)Gfu{v) and Y q {v)G fu(u) , respectively, being X q {y) and Y q (u) 
the Fourier transforms of the quadrature signals x q (t) and y q (t), and Gfu(v) = l/(l+i27ri/r) 
the (complex) gain of the filter (see figure 1). 

If the Q factor is low -as in air and in liquids- the signals (x(t)) T and (y(t)) can be directly 
used for DSFM. In fact, for low Q factors and low tip-sample interaction the frequency 
shift Aui n t induced by tip-sample interaction is smaller than the width of the resonance 
(Auint < vq/Q). Assuming the validity of the harmonic approximation for the dynamics of 
the cantilever, the signal {x(t)) T is then proportional to the frequency shift, and the signal 
(y(t)) is proportional to the oscillation amplitude, which in air and liquids is generally used 
as control parameter. For high Q factors, however, the width of the resonance is smaller 
than the frequency shifts induced by tip-sample interaction. Moreover, high Q-factors imply 
that the oscillation amplitude needs a long settling time (of the order of Q/vq) to reach its 
steady state value [17] . In this case it is necessary to track the resonance frequency using 
Phase Locked Loop techniques [17] [T8| |2T] . This is implemented with a Pi-controller that 
essentially adjusts the frequency of the voltage or numerically controlled oscillator (VCO 
or NCO, see fig. 1) so that the (x(t)) T component vanishes; the phase of the oscillation is 
then kept at —n/2 and the system is always at resonance. The output of the Pi-controller 
is then directly proportional to the frequency shift, and this is the signal used in typical 
DSFM applications. The Pi-controller represents, from an electronic point of view, a filter 
with a well defined bandwidth and gain. The time constant of the filters following the 
multiplication stages and the bandwidth and gain of the Pi-controller have to be chosen so 
that the overall closed loop system is stable [TS]. Since the precise set up of the Pi-controller 
does not modify the essential physics to be discussed here we will assume - to simplify the 
discussion - an ideal controller that instantaneously transmits variations of (x(t)) T to the 
VCO. The time response of the overall system is then determined by the time constant of 
the filters after the multiplication stage, which induce a delay of order r. 
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III. THERMAL FLUCTUATION OF THE CANTILEVER POSITION 



In addition to the "coherent" signal a ex (t) induced by the external excitation of the 
cantilever, in the present context also the "incoherent" contribution due to thermal noise, 
a t h(t), is important [18J. The total motion of the cantilever is thus a(t) = a ex (t) + a t h(t). 
When this motion is transduced into an electrical signal, the position detector will add 
some instrumental noise n(t), the total signal entering the DSFM detection unit is then 
UDSFAiit) = e (a ex (t) + a t h(t)) + n(t), where e is a constant that describes the sensitivity 
(unit:V/nm) of the photodetector electronics. The key question is now: how does this signal 
pass the DSFM detection unit and what is the final noise of the frequency output? 

The Equipartition Theorem discussed above (relation [T]) relates the total thermal dis- 
placement to the force constant of the cantilever and the temperature of the system. How- 
ever, since the DSFM - detection system performs non-trivial processing of the input signal, 
the precise spectrum of the thermal noise has to be taken into account in order to cal- 
culate the noise of the frequency output and thus the precise limit of DSFM. A detailed 
discussion of thermal noise in SFM set-ups is out of the scope of the present work (see, for 
example [91 flOl fT8l ) . In a simple picture, the noise density of the cantilever motion can be 
obtained by observing that the Equipartition Theorem (equation [TJ relates the mean energy 
of the cantilever with the thermal energy kT. Assuming that the effective thermal noise 
"force" driving the cantilever has a constant spectral noise density a t h (unit: nm/\/ Hz), 



this coupling strength and the mechanical gain G (u) of the cantilever (relation 17) define 
the spectral response of the thermal spectrum ath {v) of the cantilever: a t h (v) = ath G (v). 
The Equipartition Theorem (relation [T]) then implies: 

-kl = / av -a., (v) = - a. h / dv 



2"" Jo 2 thK 2 th J (1 - (y/ PQ f )2 + {{v/v Q )lQy 
The second integral gives v^Q tt/2, therefore the thermal noise density ath is 



, 2 kT 

a th =\ -fz 6 

The thermal noise density defines the coupling of the thermal bath into the SFM system. 
It depends on the quality factor and thus on the dissipation properties of the system. As 
the quality factor increases, thermal fluctuations induce less "thermal force" on the system, 
that is, coupling of the thermal bath and the tip-sample system becomes weaker. 



The thermal fluctuation of the cantilever position measured experimentally with a DSFM 
detection unit is calculated from the noise density in analogy to relation |5j 



"V c +bw /2 

Aa t h(v c ,bw) = a t h\l I dv \G(u)\ 2 (7) 

<v c -bw/2 

For low frequencies, that is, for frequencies well below the resonance frequency G (u) = a t h 
and the total thermal fluctuation within a bandwidth bw is Aa t h (V c > bw) ~ a t hVbw. At the 
resonance frequency G (i/) = a t hQ and the total thermal fluctuation is 



2 kT Q 

Aa th (z/ , bw) ~ athQvbw = J Vbw (8) 

V 7T c i/ 

as is well known from the literature. While the first relation is usually correct for static 
SFM applications, the second requires the condition that the measurement bandwidth is 
smaller than the width of the resonance curve: bw ~ 1/r <C Vo/Q- We note that in general 
this assumption may not be correct. In fact, for typical imaging applications (scan speed of 
1 second/line, 250 points/line) a minimum bandwidth of 250Hz is required, which implies 
Q <C 4000 for typical cantilevers (uq ~ 100 kHz) in order to fulfill the "small bandwidth" 
approximation. However, in UHV much higher Q values have been reported[25]. Note that 
the relation Aa t h (z/ ) = a t hQy/bw overestimates the thermal noise for large bandwidth, 
in fact for bw > (n/2)(u /Q) the total thermal noise according to relation [8] would be 
larger than ^kT/c, which is non-physical since it is in contradiction with the Equipartition 
Theorem (rel. [T]). As discussed below, for a sufficiently large bandwidth the total thermal 
noise is Aath = \/kT/c and is independent of the measuring bandwidth. 

A simple illustrative approximation for the spectral noise density of the tip-sample system 

is 



Q+l 
~Q 



a t h for < v < vq ( 1 



IQ 



and i/ ( 1 + 4^2 J < v < u ( 1 + ^ 



{ Q a th for i/o (l - ^) < v < u (l + ^) 
for v > i/ (l + 

This approximation shown in figure 2 satisfies as the correct noise density the rela- 
tions a„ 



tappr(^o) = Q ttappr(O), width of the peak of the order Uo/Q, total noise Aa = 
(fdis alppriv)) 1 ^ 2 = y/kT/c and, for large Q, a(0) = a t h- The relation of noise in the 
"peak", to that in the "flat" part is Aa peak / Aa fi at = Q. For high Q factors most of the 
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noise is in the resonance peak of the curve, and Aa peak « Aa to(fli = -\/kT/c. Therefore 
DSFM does not reduce but rather increase the total thermal noise as compared to static 
SFM. However, the signal to noise ratio is not modified: at low frequencies (static SFM) a 
driving force fo will induce the motion ao = /o/c, and the signal to noise ratio is ao/(a t h Vbw), 
while if this force is applied at resonance it will induce a response Qao and, for sufficiently 
small bandwidth, the thermal noise is QathVbw. Noise and signal are therefore amplified 
equally. Correspondingly, DSFM increases sensitivity by a factor of Q as compared to static 
SFM, but the (theoretical) signal to noise ratio is unchanged. Finally, we note that in nor- 
mal applications the thermal fluctuations of the x and y component are uncorrelated, and 
that both should have the same amount of fluctuation (see, however, [26]). Then, since 
(a?/»(*)>T = (\xth(t)+iyth(t)\ 2 ) T = (x 2 th (t)+yl(t)) T = (x] h {t)) T + (y 2 h (t)) T , it follows that 
(a&(*)>r = (vUt))r = (<&(*)> T A therefore 

Ax th (z/ c , bw) = Ay th (u c , bw) = Aa th (i/ c , bw) / V2 (9) 

and for the total thermal fluctuation of the x and y components: Ax t h = Ay t h = ath/ V% = 
VkT/ (2c). 



IV. FREQUENCY RESPONSE OF DSFM DETECTION SCHEMES TO THER- 
MAL FLUCTUATIONS 



To calculate the frequency noise, the relation 

dv 



Auth 



dip 



Aip 



tit 



dip 
dv 



-i 



Aip 



th 



(10) 



will be used. With cV(z/ )/cV = —2Q/uq (see relation 20) the only unknown quantity is the 
phase noise Aip t h- The phase, as defined by relation (20) is cp (t) = — 7r/2+tan _1 (x (t) jy (t)). 



We will assume that DSFM is operated in the Phase Looked Loop mode and is thus always 
at resonance, then (x(t)) r = (a th (t) /V%) T = and (y(t)} T = (Qa + a th (t) /V%) T = a os , 
with a os = Qao the oscillation amplitude at resonance. The correct calculation of the 
phase noise is non-trivial, since the phase is a non-linear function of the two variables x(t) 
and y(t), which is non- regular at the origin and thus the common rules for noise/error 
propagation have to be applied with care. The correct calculation based on statistical me- 
chanics is presented in appendix B. Here we will assume that a finite oscillation is applied 



to the cantilever in order to prevent the system to be near the origin of the {x(t),y(t)} 
phase space. Then, the mean phase is (<p(t)) = — 7r/2[27] and the fluctuation of the 
phase is Aip th = J (f 2 (t)) T - ((p(t)) 2 T = ^((tan^ 1 ^ (t) Jy (t))) 2 ) T . Using the relation 
(/ 2 (z + z)) = f 2 (z ) + \f (z )\ 2 (z 2 ) we have, with / (z) = tan^ 1 (z), z{t) = x (t) /y (t) 
and z = (x(t)/y(t)) = 0, 



th 



\ 



1 



'x(ty 



(11) 



i + (x(t)/ y (t))l \\y(t)J I _ V 

a os + a th/2) and relation|9jwe find for the phase noise measured 



And finally, with (y 2 (t)) T 
around a center frequency v c with a bandwidth bw = 1/r: 



th 



Aa 2 th (is c ,bw)/2 



(12) 



a 2 os + Aa 2 h (u c ,bw) /2 

For large oscillation amplitude, the phase noise is therefore Aa^/(v2a os ), which can be 
interpreted as a variation of the phase due to thermal fluctuation of magnitude a t h/y/2 of 
the x-component when the ^-component has an oscillation amplitude a os . For very small 
oscillation amplitude, relation 12 would give A(p t h = 1- However, as discussed above, at 
the origin of the {x(t),y(t)} phase space the phase is mathematically not well defined, 



relation (11) cannot be used, and relation 12 is not accurate. On physical arguments one 



would expect a uniform distribution of phase, that is, a (normalized) probability distribution 
p(ip) = 1/ (27r), which has a mean deviation Aip = ir/yS. As shown in appendix B this is 
correct in the limit of vanishing oscillation (see inset of figure 3). The correct relation for 
the phase noise has no simple functional relation with the oscillation amplitude, therefore 
we propose 



Aip 



th 



"th 



2a 2 os + 3Aa 2 h (u c , bw) /tt 2 
as approximation for the correct phase fluctuation which has the correct large and low 
oscillation behavior. Figure 3 shows the known large oscillation behavior for phase noise, 
the correct relation calculated in the appendix as well as the approximations according to 



relations 11 and relation 13 with the correct small and large oscillation limits. 



With relation 10 we finally obtain for the total frequency noise: 

Av t h 



Aa 2 h (i/ c , bw) 



2Q V 2a 2 os + 3Aa 2 th (u c , bw) /tt 2 



(14) 
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In order to discuss this relation, and to compare with the results known from the literature, 
we will consider the different approximations for large and low oscillation amplitudes, as well 
as for small and large bandwidth. For (very) low bandwidth, and large oscillation amplitude 
we obtain, using relation [8] 

U ll kT Q r— 

Auth = w\ 2 vbw 

2Q y it ca z ex Vq 

which is similar to the result reported in the literature [T7 ] [TS], |2"U] (relation [2]). As discussed 
in the previous section, this (very) "low bandwidth approximation" is usually not valid. 
On the contrary, we believe that in most applications the bandwidth is larger than the 
width of the resonance curve. Then, as long as instrumental noise is negligible, the correct 
approximation would be a "large bandwidth" approximation where all the thermal noise is 
"seen" by the DSFM - detection system. In this case the corresponding (total) frequency 
noise is 

A ^ = i\/^riyi (15) 

which is Au th ~ (vq/Q) a i/i/(2v / 2«ex) for large amplitude and for tc/ (2v3) Vq/Q — 0.9 u /Q 
for low amplitude. The characteristic frequency determining the thermal frequency noise is 
therefore the width vq/Q of the resonance curve, and for low oscillation amplitudes (a os <C 
ath) the thermal frequency noise is essentially given by the width of the resonance curve. 
In particular, this implies that, as demonstrated recently for spectroscopy applications [34J, 
DSFM is possible without external excitation of the tip-sample system. Moreover, we believe 
that a properly designed DFSM-electronics should be able to lock onto the thermal noise of 
the cantilever. 



V. EXPERIMENTS 

In order to confirm the validity of the relations just discussed, noise measurements have 
been made as a function of the bandwidth and the oscillation amplitude. A commercial SFM 
- system based on optical beam deflection [28] was used to measure cantilever motion and 
analysis of cantilever oscillation was performed either with the DSFM electronics [28J or with 
a digital lock-in amplifier [29]. The set-up of the SFM system and the essential features of 
the DSFM electronics are shown in figure 1 [SU] ■ A cantilever with a nominal force constant 
c ~ 0.4iV/m[3iJ was used and the tip was kept at a large (1mm) distance from the sample. 
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For this kind of cantilever, relation [T] gives a total mean (rms) fluctuation of 100pm. Figure 
4 shows the the spectral noise measurement of the cantilever movement acquired with the 
digital locking amplifier [29J. To characterize the spectral noise density and to discriminate 
thermal noise against other (technical) noise sources this data is fitted to the function 



n 



M = i + < c 

yf(i - (u/u fy + ((u/u )/qy 

The Lorentzian function is used to describe the thermal noise density of the cantilever, and 
the constant n\ ec is introduced to describe any additional (technical) noise (see also [91 135]). 
^From the fit to the experimental data, a quality factor Q = 100 ± 1, a natural frequency 



uq = 79.440 ± 0.002kHz, a thermal noise density ath=26.4±0.2 im/\/kHz and a constant 
= 17±2 im/y/kH z is found. 

The inset of figure 4 shows the total noise as a function of bandwidth, with the central 
frequency of the noise measurements at the resonance peak. In this log-log plot the square 
root dependence of the total noise on bandwidth for small bandwidth is clearly recognized 
from the slope m = 1/2. For high bandwidth, the total noise saturates. This saturation oc- 
curs for a bandwidth of the order of the width of the resonance curve Av = vq/Q ~ 0.8kHz, 



in good agreement with the discussion above (relations 12 and 15) and the data obtained 
from the spectral noise density. We note that the saturation of noise is only observed if 
other noise sources are negligible, which is clearly the case in our measurements. Then, as 
the bandwidth of noise measurement is increased only the thermal noise in the resonance 
peak is "seen" by the DSFM detection unit. If other noise sources are not negligible, then, 
as the bandwidth of the DSFM detection unit is increased (bw > v/Q), the detection unit 
will "see" this additional noise and the total noise will not saturate. Instead, it will continue 
to increase with the square root dependence known from the literature [37]. If the technical 
noise is appreciable, a total noise well above the theoretical value Aath = \JkT jc for thermal 
noise can be experimentally observed since thermal and technical noise is measured. 

Finally, figure 5 shows the noise of the frequency as a function of oscillation amplitude. 
Two different regimes are recognized: a constant regime for low oscillation amplitude where 
the total noise is independent of oscillation amplitude and a second regime where, as evi- 
denced by the slope m = — 1 in the log (noise) vs. log(a os ) plot, the noise decreases with 
the inverse of the oscillation amplitude. The transition range of this graph corresponds to 
an oscillation amplitude \JkT/c ~ 100pm, in good agreement with the value for the thermal 
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oscillation amplitude obtained from the frequency noise measurement shown in figure 4. 

VI. SUMMARY 

In the present work we have revised DSFM frequency detection and have analyzed how 
the thermal fluctuation of the cantilever is processed by a DSFM detection electronics. We 
find a general relation for the frequency noise as a function of a bandwidth and oscillation 
amplitude. This relation is correct for all possible values of parameters, while the relation 
known so far from the literature is only correct for a particular range. We find that for 
sufficiently large bandwidth -that is, small time constants of the DSFM detection electronics- 
essentially all the thermal noise of the cantilever is measured. In this case the width of the 
resonance peak is the characteristic noise of any DSFM frequency measurement for small 
oscillation amplitude, while for larger oscillation amplitude the noise decreases linearly with 
the oscillation amplitude. In the large amplitude and (very) low bandwidth limit our general 
relation converges within a constant factor, to the relation known from the literature. 

We are convinced that the results presented in this work are relevant for the precise 
determination of the ultimate limits of DSFM. In particular, our general relation shows that 
for small oscillation amplitudes the frequency noise does not diverge, but rather converges 
towards a finite value. Therefore, small oscillation DSFM might be much more competitive 
than considered up to now. DSFM without external oscillation, that is, driven by thermal 
noise, might be possible not only for the measurement of tip-sample interaction, but also 
for imaging applications. Since many high precision measurements SFM measurements -in 
particular in the field of Electrostatic and Magnetic Force Microscopy- are ultimately based 
on frequency measurements we believe that the present work will also improve understanding 
and optimization of these related SFM techniques and shed light on the ultimate limit of 
Scanning Force Microscopy in different important applications. 
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VIII. APPENDICES 
A. Appendix A 

In the harmonic approximation, the fundamental equation describing the dynamics of a 
SFM-system is that of the forced harmonic oscillator, m a(t) + 7 a(t) + c a(t) = F(t), 
where c is the force constant of the system, m its (effective) mass, 7 the constant describing 
the damping in the system and F(t) the external force driving the oscillator. With the 
definitions ujq = (c/m) 1 / 2 , Q' 1 = ^/(muo) = cool/c, and assuming a harmonic driving force 
F(t) = m cloUq cos{ui) = m cl Uq Re(e tut ) , where a is a displacement determined by the 
driving force (a = F(0)/c), this equation is transformed into 

a(t) + (uq/Q) a(t) + oola{t) = aou^ cos(wt) 

Note that, in order to avoid recurrent 2ir factors the angular frequency oj = 2nu will be 
used here instead of the frequency v as in the main text. For the steady state motion this 
equation can be solved algebraically with the classical ansatz a(t) = Re(A(u)e lu ' t ): 

Re {-u 2 A{u) e iu)t + (u /Q) iuA{u) e twt + cu 2 A(cu) e luJt } = Re(a to 2 e luJt ) 

from which the complex amplitude A(ou) is determined as 

A(u>) = ^ (16) 

1 - (uj/ujo) + i(u)/u )/Q 

A dimensionless (but complex) mechanical gain 

G(u) = ^ ( 17 ) 

1 - (uj/ujo) + i{u/u )/Q 

can be defined so that A{ui) = aoG(u). For the discussion that will follow, it is more 
convenient to describe the complex amplitude A{oj) in Cartesian coordinates: 



1 - (u/uj ) 2 



= «n., ^ , (is) 



Y(uo) = a {U 2 M/Q (19) 

(1 - (w/wb) 2 ) 2 + ((w/wo)/Q) 2 
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where X(uu) and Y(u) are the in-phase and out of phase components of the oscillation. Then 
the complex amplitude is A{oj) = X(u>) —iY(u) and the time response to the driving source 
ao^o cos(cut) is a(t) = X(uj) cos(ui) + Y{u) sin(wi) = cos (out + (p(u))) with a phase 

(fi(ou). This phase describes the delay between the driving source and the response and is 

*«) =-*--' 88)=-* / 2 +--(^r) < 2 °) 

As described in the main text, to experimentally determine state of a harmonic oscillator, 
the measured deflection a(t) is multiplied with two reference signals cos (u> r t) and sin (uu r t) 
to obtain two quadrature signals 

x q (t) = a(t) cos (ou r t) = (cos(wa£) + cos(u; s t)) + ^-^) (sin(w s t) + sin(a;Ai)) 

y q {t) = a(t) sin (u; r i) = — -— (— sin(wA^) + sin(w 2 t)) H -— (cos(u&t) _ cos(wst)) 

with wa = w - u r and ws = w + w r . With the definitions 

1 / cos(w A t) sin(w A t) \ 1 / cos(w s t) sin(w s t) 

M A (t) = - and M E (t) = - 

z \ — sin(o;At) cos(o;a£) / \ sin(cu s t) — cos(oj^t) 

the output of the two multiplication stages can be written in matrix notation as 
{x q (t),y q (t)} = (M A (t) + M s (i)) {X(u), Y(oj)}. The corresponding time evolution can 
thus be decomposed into one vector rotating clockwise with the frequency wa and another 
one rotating counter-clockwise with the frequency u^,- After the multiplication stage the 
two quadrature signals x q (t) and y q (t) are low-pass filtered. For a simple first order low pass 
the corresponding time domain signals are 



( Xq (t)) T =-f° d£x q (t-Z) 

T J — oo 



(21) 



X(u>) cos(u)At) — wa t sin(wA t) Y(ou) sin(cjA^) + wa t cos(cua t) 
~2 (l + w|r 2 ) + ~~2 (l + w|r 2 ) 

X(uu) cos(o;st) — us t sin(ws t) Y{u) sin(wEt) + we t cos(ws 
+ ^ (l + u;|r 2 ) + ^ (l + w|r 2 ) 

{y 9 {t)) T = \f Kv&-Z)* ,T ( 22 ) 

X(lo) — sin(wA^) — WA rc os(wA t) Y(uu) cos(wa^) — wa r sin(wA i) 
= ~2 (l + wir 2 ) + ^ (l + w A r 2 ) 

sin^t) + w s t cos(w s t) Y{u) - cos(w s t) + w s r sin(w s t) , , 
+ ~2 (l + w|r 2 ) + "1 (l+o; 2 r 2 ) 1 j 
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Again, two matrices 



Mr A (t) 



Mr s (t) 



2(l + r^|) 



2(1 + *M) 



1 



1 




cos(o;st) — w^T sin^st) sin(c<; 2 t) + w s r cos(o;£t) 
sin(a; s t) + w 2 r cos(a;£t) — cos(wst) + c^r sin(a;st) 



cos(cuAt) — wa t sin(a;A^) sin(u;At) + uat cos(a;st) 
- sin^A^) — Wat cos^a^) cos(o;a^) — Wat sin(a;A^) 



) 



) 



(24) 



(25) 



can be defined. The first matrix, Mta(£), corresponds to a clockwise rotation with the 
frequency wa and a delay angle <pa = — tan (cjat) while the second matrix, Mr s (t), 
corresponds to a counter-clockwise rotation with the frequency ojy, and a delay angle 
V?s = +tan (u> s t). 

B. Appendix B 

For the calculation of the variation of the phase, we will assume that the variables x(t) 
and y(t) are Gaussian variables. At resonance, as discussed above (x(t)) T = and (y(t)) T = 
Qdexc = a s, therefore their probability distributions are described by 



These distributions are normalized and have the variance Ax = Ay = a t h/y2 = 




\JkT/ (2c). To calculate the distribution of the phase, first its probability function has 
to be calculated according to the general relation (see, for example, [38J), 




which in our case leads to 
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\- fldxdy e -^I^H? e -{{y~a os )la th f 5 ^ _ tan -i ( x / y ^ 
a thJ J 

V ff ddrdr e -irsin{V)/a th f e -((rco S (V)-a os )/a th f _ ^-1 / ^ ^ , QQS /^X 

<*thJJ 

oo 

j j dfirdr e -( r2 - 2rcos W a o S +al s )/a^ h ) _ ^\ _ _ c -a? na sm 2 (ip)la th j r( ^ r e -(r-cos{ip)a os ) 2 j a\ h ) 

a thJJ ' " J 



/ oo oo 
1 e -al s siv?( V )/al h J /" rfr ( r _ cos (y,) flos ) e -(r-cos( V )a M )Va? h ) + CQS ^) ^ /" rfr g^r-co^H*) 2 /^) 



/7i 



1 -aLsin 2 M/aj / Q tfe - cos 2 ( y )^,/aj) V 



7ra 



/7i 



+ cos (</?) a os —a th ( 1 + Erf 



a>th 



COS (<£>) 



— ( e - a °° /a *» + Vtt — e - a - sin2(v)/a « cos M f 1 + Erf 
2vr V a>th V 



cos (93) 



(26) 



where Erf[x] = Ij^K^dx e~ x2 is the normalized Error Function (Erf [00] = 1). This prob- 
ability distribution is plotted in the inset of figure 3 for the range of oscillation amplitudes 
a-os/a-th = — 2. For large oscillation amplitude a os » a t h the first term can be neglected, in 
addition we can assume Erf [...] ~ 1, and only very small angles contribute to the probability 
amplitude (sin {tp) ~ <p; cos {<p) ~ 1) then, with (p th = a t h/a os , we find 



1 



nipth 



which is a normalized Gaussian probability distribution of the angle p with variance Ap 



pth/y/2 = ath/ (y/Za-os), in agreement with the high excitation limit of eq. 12. We note that 



the ratio p t h = ^th/^os can be interpreted as the fluctuation of the phase due to thermal 
variation a t h/ \/2 of the x-component when the y-component of the oscillation is fixed at a os 
(for a os » a t h thermal fluctuation of the y-component essentially gives no contribution to 



phase noise). For small oscillation amplitudes the second term in 26 is small and we find, 
to first order in l/pth, 

If 2Jtt . ,' 
PM - 7T 1 + cos VP) 

27T V p th 

The angle probability distribution therefore becomes non-Gaussian and ultimately uniform 
(see inset of figure 3), as is expected for vanishing oscillation amplitude. The mean value 
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of this (normalised) angle distribution is <p = 0, and its square deviation (Ay?) 2 = 7r 2 /3 — 
2y/nip th . The correct phase error Aip (a os /a t h) calculated from the probability distribution 



26 is plotted in figure 3, together with the different approximations discussed in this work. 
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IX. FIGURES 




Figure 1: Schematic description of a typical lock-in type DSFM detection unit. The signal 
to be analyzed by the DSFM detection is assumed to be centered around some frequency 
u Q . It enters the detection unit at the input "m", is amplified and usually high-pass 
filtered (for simplicity the corresponding components are not shown) before being 
multiplied with two reference signals in quadrature at a frequency v re f, shifting the signal 
to the frequencies vq — v re f and v$ + v re f- The resulting signals are then low-pass filtered 
to remove the higher frequency component (uq + u re f), resulting in two averaged signals 
(x(t)) T and (y(t)) . For sufficiently small interaction (x(t)) T is proportional to the 
frequency shift and can be used to re-adjust the driving frequency of the VCO (or NCO) 
by means of an appropriate feedback loop (Pi-controller). The output of the Pi-controller 
used to adjust the excitation frequency is then proportional to the frequency shift 8v(t). 
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Figure 2: Simple approximation of the noise density (black) for a thermally excited 
cantilever as discussed in the main text together with the correct noise density (red). For 
large quality factors, most of the noise is within the main peak at the resonance frequency. 
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Figure 3: Main graph: Thermal noise error of the phase as a function of the (relative) 
oscillation amplitude a os /a t h- The black, solid, thin line corresponds to the relation 
obtained from the relation known in the literature, which diverges for small oscillation 
amplitude. The black, dotted line corresponds to the relation Atp t h = a/ ((tan _1 (x/?/)) 2 ) 



(11), which is not correct at the singular point {x,y} = {0,0}. The red, thick, solid line 



shows the correct relation calculated from the probability distribution 26 discussed in the 
appendix B. Finally, the red, thin, dotted line corresponds to the approximation 



Aipth = \J ctffj {1a? os + 3a 2 h / it 2 ) (relation 13), which has the correct low and large amplitude 
limits. 

Inset: Probability distributions p v ((f) for different (relative) oscillation amplitude a os /a t h- 
The probability distributions have been calculated for the range of oscillation amplitudes 
a s/ a th = — 2. The probability distribution p^i^p) for a os /a th = is flat while that for 
a os/ a th — 2 is essentially Gaussian and has the highest peak at (p = 0. 
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Frequency [kHz] 



Figure 4: Main graph: Spectral noise density of a 0.4N/m cantilever measured with a 
digital lock- in amplifier. For this noise measurement, no external excitation was applied to 
the cantilever and the motion of the cantilever was measured using the beam-deflection 
technique. The larger (red) points correspond to experimental noise data, the solid line to 
a fit assuming a constant offset and a Lorenz function (see main text) and the smaller 
(pink) points show the error between this fit and the measured data points. Inset: Log-Log 
plot of the total noise as a function of bandwidth for a noise measurement centered at the 
peak of the main noise curve. For small bandwidth, the frequency noise shows the typical 
l/a os behavior (slope -1 in the Log-Log plot). However, for high bandwidth (bw > fo/Q, 
with Q quality factor) the total noise saturates. 
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Figure 5: Frequency noise of the DSFM detection electronics measured as a function of 
oscillation amplitude for two different bandwidths (50 Hz and 100 Hz). For this 
measurement, the same cantilever as that used for the previous experiment was utilized 
(force constant of 0.4N/m). The cantilever was excited by the DSFM electronics with the 
phase-locked loop enabled, and the frequency output Su(t) was fed into a digital lock-in 
amplifier in order to determine the total noise of the frequency measurement of the DSFM 
detection unit. For small oscillation amplitude of the cantilever (a osci < a t h, see main text), 
the frequency noise is independent of oscillation amplitude and for large amplitude the 
noise decreases linearly (slope 1 in the log — log plot). 
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